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Abstract. We compute the 0(1/Nj) corrections to the flavour non-singlet quark bilinear 
currents in QCD in arbitrary spacetime dimensions. Hence, the anomalous dimension of the 
QED current ipa^ u ip is deduced at four loops in the MS scheme up to one unknown parameter. 
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Recently, the large Nf method of computing the renormalization group functions of a quan- 
tum field theory has been extended to the determination of information in quantum chromo- 
dynamics, (QCD), at a new order of expansion, 0(1/ Nf), Q. In [|, § the critical exponents 
corresponding to the anomalous dimensions of the quark field and the quark mass were computed 
in d-dimensions by exploiting the structure of the field theory at the d-dimensional infrared sta- 
ble fixed point of the /3-function as well as using the relation between QCD and the reduced 
model known as the non-abelian Thirring model, (NATM), ||. It was demonstrated in H, |2| that 
the gluon field strength operator became irrelevant at the fixed point allowing one to perform 
calculations with only the quark gluon vertex and the ghost vertex. This important observa- 
tion, ||, paved the way for the new 0{1/Nf) calculations in QCD. Whilst the quark anomalous 
dimension is a gauge dependent quantity, its evaluation in d-dimensions in the Landau gauge 
is necessary for other 0(1/ Nf) computations, since like in explicit perturbative calculations the 
wave function renormalization needs to be performed first. However, in additionally providing 
the d-dimensional value of the quark mass anomalous dimension, which is a gauge independent 
quantity, it was possible not only to verify the correctness of the recent four loop MS perturba- 
tive quark mass dimension, Q |j, at 0(1/Nf) but also to determine several of the coefficients 
which will appear in the five loop MS anomalous dimension. Indeed given the huge complexity 
of performing the very high order computations in QCD, insight into the large order structure 
of the renormalization group functions can currently, we believe, only be gained by approaches 
such as the large Nf expansion coupled with other techniques. For instance, a recent study by 
Broadhurst, ||, of QED in the quenched approximation involving the Schwinger Dyson equation 
has allowed several new coefficients to be determined in the three and four loop MS anomalous 
dimensions of various operators. By quenched we mean that part of the renormalization group 
functions which does not depend on Nf and hence arises from Feynman diagrams with no elec- 
tron (or quark in the case of QCD) bubbles. Clearly given the amount of resources required 
for large loop calculations it would seem that the construction of the high order terms could 
be achieved by the marriage of the large Nf approach which performs the electron or quark 
loop bubble sums with a method which focuses on the harder Nf = graphs which have no 
subdiagrams. Therefore, it is the purpose of this article to extend the 0(1/ Nf) calculation of 
of the anomalous dimension of the (flavour non-singlet) quark currents ^Aip, where A rep- 
resents a set of 7-matrices, to 0(1/Nf). Since the case .4 = 1 corresponds to the quark mass 
operator, our calculation builds on the 0(1 /Nf) d-dimensional computations of fll], |2|. As a 
consequence of the result we present here, it will be possible, for instance, to deduce three of the 
four MS coefficients in the polynomial in Nf of the anomalous dimension of the QED operator 
tjja^ip, where a^ u = 5 [7^, 7^], with the quenched result of Q. Moreover, since the three loop 
MS anomalous dimension of the general set of operators ipAtp has been provided in explicit 
perturbation theory in QCD in Q, there is a non-trivial perturbative check on our final result. 
A final motivation for considering the renormalization of these operators rests in the fact that 
they underly the running of various quantities in heavy quark effective field theory, 0, and it 
is therefore important to have information on their anomalous dimensions. 

We will determine the anomalous dimension of the set of operators 

(n) = (1) 

where p^ 1 ---^™ i s defined by 

(n) J 

pAtl.-Mn = 7 tW _ _ _ ^Mn] (2) 

and is totally antisymmetric in the Lorentz indices for n > 1. This particular combination of 
7-matrices is chosen since there will be no mixing under renormalization and also because the 
set {r(„)} corresponds to a complete 7-matrix basis in (^-dimensions, ||. The properties of the 
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matrices I\ n ) are well established and given in flu], 11, 12]. To compute the large Nf contribution 



to the operator anomalous dimension one follows the standard critical point approach of |lj: 
adopted for QCD, ]i~4] , g]. This involves using the NATM where there is no triple or quartic gluon 
self interactions and the quark and gluon propagators are replaced by their leading d-dimensional 
critical point forms, 
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where A and B are the momentum independent amplitudes, a and (3 are the respective dimen- 
sions of the quark and gluon fields and b is the covariant gauge fixing parameter. The field 
dimensions are related to their anomalous dimensions via 
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where rj is the quark anomalous dimension critical exponent, \ ls the anomalous dimension of 
the quark gluon vertex and the spacetime dimension d is related to /i by d 



For the 

Feynman diagrams which arise at the order we consider none involve ghost fields which is why 
we have omitted its critical propagator. Ghost contributions are present indirectly in the final 
critical exponent for the operators 0( n ) through the contribution from rj and the combination 
z = A 2 B which arises from each vertex. Then the operators 0( n ) are substituted into the two 
point Green's function (V ; (p)C(n)V ; ( — p)) an d the scaling behaviour of the appropriate diagrams 
determined using (||). In particular a regularization with respect to A is introduced by shifting 
the anomalous dimension of the vertex, (3 — > f3 — A, and the poles with respect to A are removed 
in a standard renormalization by the appropriate renormalization constants. If p is the external 
momentum then the resulting finite Green's function will involve a lnp 2 term. For that Feynman 
diagram its coefficient contributes to the overall anomalous dimension critical exponent of the 
operator, |15[| . By summing all such contributions at each order in 1/Nf one arrives at the full 
anomalous dimension. As we are interested in the 0(1/ Nf) calculation there are only seven 



Feynman graphs to consider which are illustrated in figures 1 and 2 of 1 16] as well as the three 
loop graphs of figure 2 of Q where the operator insertion is on a fermionic line which joins to 
an external vertex. 

Although it was argued in [Q] that when computing in the critical point large Nf approach, 
the Landau gauge was the only sensible choice of gauge, we have chosen to calculate the operator 
dimension in the Feynman gauge. The main reason for this is that since the operators 0( n ) are 
gauge independent, their anomalous dimensions are therefore independent of b and hence can 
be determined in any gauge. We therefore set 6 = 1 in @ to reduce the number of terms in the 
Feynman integrals and thereby minimise the amount of computation to be performed. Further, 
we will require rj and z at 0(1 /Nf) in this gauge. Since we are following the original skeleton 
Schwinger Dyson approach of |l5|], we record their respective values by first setting 
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Then for arbitrary 6, Ji~4| , g], 

[(2 M -1)( M - 2) +bn]C 2 {R)ito 



r<> + l)r] 



+ 



(2h-1)(jm-2)T(R) ' ~ ! 2(2 / j-l)( M -2)T(i?) 

+ nb - 5/x + 2)Gu - 1) [£ 2 0u) + $(/i)] C 2 (G) 

(8/x 5 - 92// + 270n 3 - 301/x 2 + 124/i - 12)(2/i 2 + fxb - 5/x + 2)j>(n)C 2 (G) 

2(2 fj, - l)(2/i - 3)(/x - 2) 



+ 3/x(/i - l)\pLbC 2 {G) + 2(2/i 2 + /ib - 5/x + 2)C 2 (i?)]e(/x) 

- [(32// 7 & - 96// 7 + 8/x 6 6 2 - 224// 6 6 + 912// - 4// 5 6 3 - 84// 5 6 2 + 704//6 - 3360// 5 

1124/i 4 6 + 6240// 4 - 19/i 3 6 3 - 387// 3 6 2 + 846/i 3 6 - 6292// 3 



+ 16/i 4 6 3 + 278// 4 6 2 
+ 6/t 2 6 3 + 230/x 2 6 2 - 



222// 2 6 + 3416/i 2 



48//6 2 - 4//6 - 908// + 88)C 2 (G)// 



8(4// 5 + 4// 4 6 - 32/ - 13// 3 6 + 75// + 8// 2 6 



70/i 2 - 2//6 



+ 32// - 6)(2/i - 1)(2// - 3)(// - 2)C 2 (R)] 
/[4(2//-l)(2//-3)(//-l)(//-2)//] 



C 2 (i?)r? 2 



2(2/i - 1) 2 (// - 2) 2 T 2 (R) 



(6) 



where the second result derives from the 0(1/N 2 ) solution of the Schwinger Dyson 2-point 
equations for non-zero b, and 



Z-2 



3^-l)T(ti + l)C 2 (R)rj 
2(2//-l) 2 (//-2) 2 T 2 (i?) 

M 2 ( M -l)r(/x)(7 2 (G)r 7 2 



+ 



e( M ) - 

39(//) 



1 



4(2//-l) 2 (//-2) 2 T 2 (.R) 

(8// 5 - 92// 4 + 270// 3 



(//-l) 2 
- * 2 (m) 



6 2 



4(//-l) 2 



+ 



301/i 2 + 124// - 12)*(//) (/x 2 + 2/x - 2)6 



2/i(/i - 1)(2// - l)(2/i - 3)0 - 2) 
(16// 7 - 120// 6 + 420/i 5 - 776// 4 + 742// 3 



+ 



2(/x-l)V 
349/i 2 + 84// - 12) 



2(2// - l)(2/x - 3)(// - l) 2 (/i - 2)/x 2 



(7) 



where 
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and ip(x) is the derivative of the logarithm of the Euler T-function. 

Since the evaluation of the 0(1/ N 2 ) Feynman diagrams follows the standard methods of 
0, 11, 14], we quote the final value of the critical exponent for the operator C?( n ) for arbitrary 
n. We find 
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in the notation of (@). The value for r/j n ^ was given in 
in our leading order calculation. 



and we have correctly reproduced it 
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There are various checks on the result, @. First, since the operator -ip^ip is a conserved 
physical current, its anomalous dimension must vanish at all orders in perturbation theory. The 
overall factor of (n — 1) which naturally emerges in (||) ensures this. Second, the result for the 
quark mass anomalous dimension at 0(l/iVj), jl], is recovered when n = 0. The remaining 
checks are a consequence of comparing with explicit perturbative calculations since the three 
loop MS anomalous dimension for Or n ) has recently been provided in || and is 
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where (,{z) is the Riemann zeta function and a is related to the strong coupling constant, a s , by 
a = a s /(47r). The lower order coefficients for n = 2, for example, had been established in jj], |l7]] 
whilst the quenched term at three loop in QED had been given in Q. To compare @ with ( JlOj ) 
we evaluate the latter at the critical coupling, a c , in d-dimensions and expand in powers of e 
and 1/Nf, where d = 4 - 2e, with Q 



3e 



T(R)N f 



99 



+ 



1 



T 2 {R)Nf 
237 



33 



27 



45 



— C 2 (G)e - — C 2 {R) + —C 2 (G) e 



16 



32 



tt;C2{R) + ^rC 2 (G) e 3 + — C 2 (i?) + — C 2 {G) e 



77 



53 



16 



32 



^ [(288C(3) + 2U)C 2 (R) + (480C(3) - 229)C 2 (G)] + 0(e 6 ) 



To 0(e 3 ) we find exact agreement between both expansions of the critical exponent. Moreover, 
since (10) was computed in an arbitrary covariant gauge and (|9|) is in agreement with it, we 
have justified the choice of Feynman gauge in our calculations. 

Having established the correctness of @ with all regions of overlap, we can now determine 
new information on the higher order terms of (|To|). Writing the anomalous dimension in terms 
of its 0(1/Nf) part as 

7(n)(a) = - (n - l)(n - 3)C 2 (R)a 

+ U(n - 15)T(R)N f + (18n 3 - 126n 2 + 163n + 291)C 2 (G) 
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and expanding (0) to 0(e 5 ) using (|TT]) , we find the new MS coefficients 
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c 52 = (n - 1) [96(72n 3 - 504n 2 + 17n + 2013)C(3) - 3456(lln - 51)C(4) 

+ 110592(n - 3)C(5) - (7128n 3 - 37416n 2 + 64723n - 56001)]/729 . (13) 

As part of our motivation was to make use of other techniques aside from explicit perturba- 
tion theory to deduce higher order terms of the renormalization group functions, we can now 
examine the operator 0/ 2 ) in QED at four loops. The anomalous dimension in the quenched 
approximation is available from Q and restricting fll3| ) to QED, we find that 
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where only one unknown coefficient, C43, remains to be determined. At present there is no next 
to quenched approximation technique which would allow for it to be deduced aside, of course, 
from an explicit calculation at four loops where one would isolate only those Feynman graphs 
which were linear in A^. 

Finally, we address the issue of the relation of the anomalous dimensions of Oua to their four 
dimensional counterparts. Due to the arbitrary dimensional nature of the large method, we 
were forced to use the infinite dimensional T( n ) basis which when restricted to four dimensions 
is related to the conventional 7-matrices. However, one has to be careful in treating 7 s which 
is a strictly four dimensional object and never arises in the generalized basis. For example, 
setting n = 0, 1 and 2 in (^) we recover the correct perturbative structure of the respective 
flavour non-singlet currents. For n = 3 and 4 the anomalous dimensions r]^ and r]^ do not 
correspond to those of the axial vector and pseudoscalar currents, [18]. The reason for this is 



that the operator with T^), for example, does not exactly project onto the operator with 7 s in 
four dimensions and in some sense incorrectly retains an anticommuting 7 s in d-dimensions. In 
the large A^ approach the method to properly account for the presence of 7 s has been developed 
in |lj| based on earlier observations of [20]. Like perturbation theory, []l8||, one has to add an 
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additional critical exponent which would correspond to a finite renormalization constant in the 



perturbative approach. However, the result of [19 1 only concerned 0(1/Nf) calculations and 
therefore we extend the method here to the next order since it will also be important in future 
large Nf calculations of other operators which contain a 7 s such as four quark operators. As in 



[19] we define a critical exponent corresponding to the finite renormalization by 

V &n = r/ 4 ~ n ) - r/ n ) . (15) 



This choice is motivated by [15, M in that since the currents are non-singlet, and hence anomaly 
free, this condition will restore the four dimensional properties of 7 s in the renormalization of 
the operators. Hence, from @ we find, 
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which is only relevant when n = 3 or 4. Adding ( |i~6| ) to @ for these values one correctly recovers 
the critical exponents for the axial vector and pseudoscalar currents which agree with the known 
four dimensional perturbative results, [18]. 

We conclude with various remarks. First, we have demonstrated that it is possible to marry 
results from several techniques aside from perturbation theory to determine the structure of a 
renormalization group function at four loops in MS, up to one unknown parameter. Second, 
we focused on the flavour non-singlet currents. For even n the results are also valid for the 
flavour singlet currents since the extra graphs which arise when the operator 0( n \ is inserted in 
a closed quark loop trivially vanish as they involve a trace over an odd number of 7-matrices. 
The study of the flavour singlet currents at 0{1/Nj) for odd n is more involved than the non- 
singlet computation. Due to the nature of the 1/Nf expansion the largest loop order one would 
have to consider is five loops. However, similar diagrams arise in the computation of the quark 
gluon vertex anomalous dimension at 0(1/ ) and since that more fundamental calculation has 
not yet been performed, the formalism and integration techniques necessary to handle the five 
and lower loop diagrams in the critical point approach using the propagators (||) are not yet 
available. Therefore, to extend @ to the singlet case, prior to treating the problem of correctly 
accounting for the axial vector anomaly, which is an issue in its own right, a substantial amount 
of additional calculations would be required which is beyond the scope of the present article. 
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